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Abstract 

Coherent resonance is the effect of resonant excitation of nonhnear coherent modes 
in trapped Bose condensates. This novel effect is shown to be feasible for Bose- 
condensed trapped gases. Conditions for realizing this effect are derived. A method 
of stabilizing Bose condensates with attractive interactions is advanced. The origin of 
dynamic critical phenomena is elucidated. Interference effects are studied. The exis- 
tence of atomic squeezing and multiparticle coherent entanglement is demonstrated. 
Coherent resonance is a generalization of atomic resonance, involving internal states of 
an individual atom, to collective states of a multiparticle system. 
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1 Introduction 



Resonant interaction of an alternating electromagnetic field with an atom makes it possible to 
select a pair of atomic levels and to treat the atom as an effectively two- level system [1]. This 
resonant interaction is in the basis of one of the main directions of optics. The possibihty 
of realizing atomic resonance is due to the existence of discrete energy levels of electrons 
in an atom, the energy spectrum being described by the Schrodingcr equation. The Bose- 
condensed atoms are also described by the Schrodinger equation, though nonlinear, which 
is often termed the Gross-Pitaevskii equation [2-4] . Trapped atoms, because of a confining 
potential, also possess a discrete spectrum [5]. A specific feature is that this is a spectrum 
of collective nonlinear states of an ensemble of coherent atoms. The collective coherent and 
nonlinear nature of these states is what makes them principally different from single-particle 
linear states of an individual atom. To emphasize this difference, the collective states of Bose- 
condensed atomic gases are called nonlinear coherent modes [5] . Since the spectrum of these 
modes, pertaining to trapped atoms, is discrete, it has been suggested [5] that a resonant 
coupling of two such modes is feasible by means of an alternating field with a frequency tuned 
to the transition frequency between the chosen levels. This type of resonance for exciting 
nonlinear coherent modes can be named coherent resonance. This phenomenon, if realised, 
could open a whole new field of possible applications, analogous to those accomplished with 
atomic resonance. 

In the present communication, we investigate the conditions required for realizing the 
coherent resonance. We propose a new method for stabilizing Bose-condensed gas with 
attractive interactions by transferring atoms to an excited coherent mode. We study dy- 
namic critical phenomena, discovered numerically [6-8], and give their explanation. We also 
describe some new properties of Bose condensate subject to coherent resonance, such as 
interference effects, atomic squeezing, and multiparticle coherent entanglement. 

2 Coherent Resonance 

We consider a system of trapped Bose-condensed atoms at low temperature, when the Bose 
gas in a coherent state. The coherent-field wave function satisfies the Gross-Pitaevskii equa- 
tion 



in which the coherent field is normalized to unity, ||<^|p = 1. The nonlinear Hamiltonian is 




where U (r) is a trapping potential; mo is atomic mass; N is the number of atoms; Ag = 
Anh'^as/mo, with being a scattering length. The alternating field is 




(1) 







(3) 
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Note that Eq. (1) is an exact equation for the coherent field [9], and v^(r,t) does not need 
to be interpreted as an average of a field operator, for which case Eq. (1) would be only a 
mean-field approximation. 

The nonlinear coherent modes arc the stationary solutions to the eigenproblem 



(4) 



Since the Hamiltonian (2) is nonlinear, the modes are not necessary orthogonal to each other, 
so that the scalar product 

{<Pm:<Pn) = j <^^(r) <^n(r) dr 

is not compulsory a Kroneker delta. But each mode can always be normalized, ||<^n|P = 1- 
Some properties of the nonhnear modes have been discussed [5-8,10-13] and dipole mode 

has been observed experimentally [14]. 

Let us select two spectrum levels, with the energies Ei and E2, such that Ei < E2. The 
related transition frequency is 

UJ21 = J (^2 - El) . (5) 

The frequency of the alternating field (3) is tuned close to this transition frequency (5), 
which implies the resonance condition 



Aa; 



•C 1 , Aa; = uj — UJ21 ■ 
The solution to the evolution equation (1) can be presented in the form 

fi^^t) = J2^n{t) <^n(r) exp ^Ent^ . 

We assume that c„(t) is slow as compared to the exponential, such that 

dcn 



(6) 



(7) 



En 



dt 



< 1 . 



(8) 



There are two transition amplitudes in the system, one is the internal transition amplitude 

N 



amn = As — J \(Pm{r)\'^ [2|(^„(r)|2 - |<^m(r)|^J dr , 
caused by atomic interactions, and another is the external transition amplitude 

Pmn^l J ¥';(r) mr) - iV2{r)] <p^{r) dr , 



(9) 



(10) 



due to the driving resonant field. Similarly to atomic resonance, in order to avoid power 
broadening, the transition aplitudes are to be smaller than the transition frequencies. 



«1, 



< 1 . 
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Under conditions (6), (8), and (11), the Gross-Pitaevskii equation (1) reduces to the set of 
equations 



^ -TT = "12 C2 Ci + - i3i2C2e 

at 2 



• I ,2 1 ^* 

I — = a2i\ci\ C2 + - Pi2Cie 



iAuit 



;i2) 



for an effective two-mode system, where |cip + |c2p = 1. This reduction is vahd provided all 
transition amplitudes, involved in Eq. (12), are small in the sense of condition (11), which 
yields the inequalities 



"12 



<^21 



«1, 



"21 



<^21 



«1, 



12 



<^21 



< 1 . 



(13) 



To analyse the above inequalities, we consider a cylindrical trap modelled by the harmonic 
trapping potential 



with the aspect ratio 

Introduce the dimensionless coupling parameter 
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(14) 



(15) 



The spectrum of nonlinear coherent modes is defined by the eigenproblem (4) . The external 
transition amplitude (10) can always be made sufficiently small by regulating the amplitude 
of the resonant field. Wc need to calculate the internal transition amplitude (9) and to 
compare it with a chosen transition frequency (5). Calculations can be accomplished by 
means of the optimized perturbation theory [15-17] (for reviews sec [18,19]). The first two 
inequalities (13) can be valid only outside the radius of convergence of the strong-coupling 
expansion. This imposes the restriction 



\9^\ < 



(16) 



nmj 



on the coupling parameter (15). Here the integral Inmj ~ (I'^nmjP, iV'nmjP), and p = 2n + 
|m| + 1, q = 2j + 1 are the combinations of the radial, n, azimuthal, m, and axial, j, quantum 
numbers. The restriction (16) can be rewritten as the limitation on the admissible number 
of particles N < Nq, with the limiting number 



[2p^ + {quYf/^ 
dGnupy/q Inmj 

For the ground-state level, this reduces to 

/7f (2 + Z/2)5/4 



(17) 



;i8) 
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Thus, the number of particles, which can be resonantly transferred to an excited coherent 
mode, is limited by a number Nq, which depends on the quantum numbers of the coupled 
modes. The number Nq also essentially depends on the trap shape through the aspect ratio 
(14). 

It turns out that the limiting numbers (17) and (18) are close to the critical values 
defining the maximal number of atoms with a negative scattering length, which can be 
condensed in a trap. Examples of atoms with attractive interactions are ^Li (see review 
[20]) and ^^Rb (see Ref. [21]). Such atoms can form condensates only if their number 
does not exceed a critical value [5,22-25]. In that case, the condensate forms a long-lived 
metastable state, with a very slow decay, caused by quantum tunneling [26-28], with the 
tunneling probability being negligible as compared to the probability of escaping from the 
trap because of depolarizing collisions [23,29]. But if the number of atoms is larger than 
critical, the condensate collapses. With a supply of atoms from an external source, say by a 
continuous loading from another trap [30], the condensate grows again and, after surpassing 
the critical number, again collapses. Thus a series of growths and collapses take place, as 
was observed experimentally [31,32] in Bose-condensed ^Li and described theoretically by 
means of the Gross-Pitaevskii equation supplemented by relaxation terms and by the related 
rate equations. Explosion of an attractive condensate of ^^Rb atoms by manipulating the 
atomic interactions with an external magnetic field near a Feshbach resonance [36] has also 
been observed [21,37]. 

To increase the number of atoms in a condensate with attractive interactions, it was 

suggested [38] to drive a quadrupole collective excitation. We may notice that transferring 
atoms to an excited nonlinear coherent mode may also stabilize such a condensate. This 
follows from the form (17), which shows that the limiting number Nq increases for higher 
excited modes as 

A^o ~ (2n + |m| + 1)=^/' , A^o ~ (2j + 1)' . 

To estimate the number of atoms for highly excited modes, one could also use an optimized 

quasiclassical approximation [39]. As Eq. (18) demonstrates, the limiting number Nq is 
larger for essentially anisotropic traps, that is, for cigar-shape {v ^ 1) and disk-shape 
{u ^ 1) traps. 

The temporal behaviour of the system under coherent resonance is described by the 
evolution equations (12). Numerical solution of these equations, for the varying amplitude 
of the resonant field, displays dynamical critical phenomena [6-8] . Here we present a general 
picture of these phenomena and elucidate their origin. 

Let us present the population amplitudes in the form 




C2 = y exp |i [n2- ^ ^) } ' (19) 

where tti = 7ri(t) and 7r2 = 7r2(t) are the phases. The variable 

s=|c2p-|cip (20) 

is the population difference. Introduce the notation for the average internal transition am- 
plitude 

ot = - {ai2 + a2i) (21) 
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and let us take into account that the external transition amplitude (10) is, in general, complex 
valued, 

Pu^Pe'\ P=\f3i2\. (22) 
Also, define the effective detuning 

S = Au; + ^ (ai2 - 021) • (23) 
And introduce the phase difference 

X = 7r2 - TTi + 7 . (24) 

Substituting Eqs. (19) into the evolution equations (12), we come to the set of equations for 
the population difference (20) and phase difference (24), 

ds ^ r- ^ . dx i5s - 

= — pv 1 — s smx , — — as ^ , cosx + . (25) 



dt ' dt y/T 



s 



Numerical solution of Eqs. (25) demonstrates the existence of dynamical critical phe- 
nomena, when a tiny variation of parameters results in sharp changes of temporal behaviour 
of s{t) and x{t). We illustrate this in the following figures, where time is measured in units 
of a~^. There are two independent dimensionless parameters 

b=^, £ = -, (26) 
a a 

which can be varied. In what follows, we vary h in the interval [—1, 1] and keep |£:| <^ 1. 
And there are also two initial conditions sq = ^(O) and xq — x(0). 

In Fig. 1, we set sq = ~1) — 0, and e = 0, while varying the parameter h corresponding 
to the amplitude of the pumping resonant field. In the range < 6 < 0.5, the population 
difference oscillates between —1 and zero and the phase difference monotonically decrease, 
as is shown in Fig. la and lb. The critical point he = 0.5 separates the regions of two 
qualitatively different temporal behaviours. After h surpasses he, the population difference 
starts oscillating between —1 and +1, while the phase difference becomes a periodic function, 
as is demonstrated in Figs. Ic to le. 

As follows from the definition of the internal transition amplitude (9), its value may be- 
come negative, either for strongly energetically separated modes or for atoms with attractive 
interactions, when Ag < 0. In such a case, the parameter (21) is negative, a < 0, and so 
is the parameter & < from Eq. (26). The variation of h in the negative region is also 
accompanied by dynamic critical phenomena, as for positive b, though the behaviour of the 
phase difference is not the same as for 6 > 0. This is illustrated by Fig. 2, where we set 
the same initial conditions sq — —1 and Xq = 0, and zero detuning. For b in the interval 
—0.5 < 6 < 0, the population difference is in the range — 1 < s < 0, as is shown in Figs. 2a 
and 2b. The critical point now is be = —0.5. After crossing the critical point, when b < be, 
the amplitude of the population difference increases by a jump, so that s{t) oscillates now 
in the diapason — 1 < s < 1. The phase difference is an oscillating function for both b > be 
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and b < be, but its shape changes drastically when crossing be, as can be seen in Figs. 2c 
and 2d. 

In this way, there are two critical lines for each given pair of initial conditions Sq and Xq. 
Changing initial conditions changes the location of the critical line on the manifold of the 
parameters b and £. Thus, in Fig. 3, we show the course of events for the initial conditions 
So = —0.8 and Xq = 0, with the detuning e ~ and for the parameter b varying in the whole 
range —l<b<l. Under the fixed detuning £ = and given initial conditions, there are 
two critical points bd = —0.8 and bc2 = 0.2. 

For a nonzero detuning, the location of the critical line changes, but the overall picture 
remains the same. A finite detuning also makes the shape of the function x{t) a little 
asymmetric, as is shown in Fig. 4. 

To understand the origin of these critical phenomena, we studied the phase portrait on 
the plane {s,x} for different parameters b and e. Examples are demonstrated in Fig. 5. We 
also accomplished the stability analysis of fixed points. This investigation clarify the origin 
of the dynamic critical phenomena. The latter are related to the existence of the saddle 
separatrix defined by the equation 



The separatrices divide the phase plane {s, x} into the regions with different dynamics. As 
is evident from Eq. (27), changing any of the parameters 6 or £ shifts the separatrices. When 
a separatrix crosses an initial point {so,Xo}, the system trajectory passes to another region 
of the phase plane, as a result of which the system dynamics changes qualitatively. The 
separatrix crossing of an initial point corresponds to the critical line 



on the manifold of the parameters b and £. It is this separatrix crossing effect that causes the 
appearance of the dynamic critical phenomena. It can be shown [7,8] that a time-averaged 
system displays critical phenomena, typical of phase transitions in stationary statistical 
systems, on the same critical line. 

There exist several more interesting effects occurring in the process of coherent resonance. 
Not to overload the present communication, we mention these effects only briefly. A more 
detailed consideration will be done in separate publications. 

2.1 Interference effects 

Coherent resonance couples two selected nonlinear modes, leaving other modes unpopulated. 
Therefore, the coherent field (7) is, effectively, a sum 



b\/l — s"^ cos X -\- £s — b . 



(27) 




(28) 



ip{v,t) = ipi{v,t)+ip2{v,t) 



(29) 



of two terms 
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Using this, one can define the interference pattern 



Pint{r, t) = p{r, t) - pi(r, t) - p2{r, t) , (30) 

in which 

p{r,t) = \^{r,t)\^ , pi{r,t) = \^i{r,t)\^ , 
and the interference current 

jmt(r, t) = j(r, t) - ji(r, t) - j2(r, t) , (31) 

where 

h h 

j(r,^) = — Im ip*{r,t)Vip{r,t) , ii{r,t) = — Im ip*{r,t)V(pi{r,t) . 
mo mo 

These quantities show that the total density of atoms in a trap is not simply a sum of 

the partial mode densities but it includes also interference fringes and that, because of the 
different mode topology, there arises a local interference current. The interference effects 
can be observed experimentally. 

2.2 Atomic squeezing 

The evolution equations for the effective two-mode system, obtained by means of the coher- 
ent resonance, can be interpreted as the equations for the statistical averages of quasispin 
operators 

N 

Sa = Sj^ , S± = Sx i iSy , 

i=l 

where a — x,y,z and Sf is an a-component of the spin-1/2 operator. For instance, the 
population difference (20) can be written as the average 

2 

N 

The squeezing factor may be defined [40] as 

^< >2 + < 5„ >2 I < S± > I 

where A'^{Sz) =< Si > — < One says that is squeezed with respect to S± if 

Qz < 1. Calculating the factor (32), we find 



Q, = . (33) 

Since < 1, one has < 1. The maximal squeezing is achieved when one of the modes 
is completely populated, i.e. s = ±1. Squeezed 5*2 means that the population difference can 
be measured with a better accuracy than the phase difference. 
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2.3 Multiparticle entanglement 



The system of interacting atoms possessing internal states may form multiparticle entangled 
states [41,42]. In our case, atoms have not internal but collective coherent modes. This 
makes the principal difference of the resonant Bose condensate having collective nonlinear 
modes from atoms with internal single-particle states. But the resonant Bose condensate 
also forms a multiparticle entangled state. This follows from the Schrodinger representation 
for the wave function of N atoms in the resonant condensate. 



N N 



*jv(ri, r2, . . . , r^, t) = Ci(t) [] Vii^i) + C2{t) ^2(^1) , (34) 



1=1 i=l 



where v^j(r) is a coherent j'-mode, with j = 1,2. The function (34) cannot be factorized 
into a product of single-particle functions, provided that ci 7^ and C2 7^ 0. The maximal 
entanglement is reached for |ci| = \c2\ — l/V^, that is, when s = 0. Since the coefficients 
ci{t) and C2(t) are defined by the evolution equations (12) or (25), from where 



,2 I |2 1 + S 



|cir = , \C2\' = , (35) 

their values can be manipulated by applying the resonant field and switching it off at the 
appropriate moment. In this way, it is possible to create any degree of entanglement and 
then one can disentangle the created entangled state. The resonant field, thus, acts as a 
quantum disentanglement eraser [43]. 



3 Conclusion 

In trapped Bose condensates, the effect of coherent resonance is feasible. This is achieved by 
applying an alternating field with the frequency tuned to the transition frequency between 
two nonlinear modes. To preserve the resonance condition requires that the number of atoms 
be limited by a maximal value. Estimates show that this value is around 10^ for a typical 
spherical trap and for alkali atoms. The limiting number is higher for cigar-shape and disk- 
shape traps, and can be about 10^ atoms. This number also increases for higher excited 
nonlinear modes. Temporal behaviour of fractional populations exhibits dynamic critical 
phenomena occurring on a critical line in the parametric manifold. The origin of these 
critical phenomena is the saddle separatrix crossing by the starting point of a trajectory. 
Actually, there exists a whole bunch of critical lines, each of which corresponds to a different 
starting point. 

As a result of coherent resonance in a system of Bose-condensed atoms, a new type 
of matter is formed, which can be called resonant Bose condensate. This is analogous to 
a resonant atom, with the principal difference that the resonant condensate is a coherent 
multiparticle system. Collective nature of the latter causes the appearance of interference 
patterns, interference current, atomic squeezing, and of multiparticle coherent entanglement. 

We have considered here a single-component Bose gas. Possible extensions of the theory 
could be to binary mixtures of Bose condensates and, generally, to multicomponent con- 
densates composed of different atomic species. Another possible generalization could be to 
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the systems of cold atoms in optical lattices [44] and to radiating atoms interacting with 
electromagnetic field [45] . We think that the effect of coherent resonance can find numerous 
applications, such as creating selected modes of atoms lasers, information processing, and 
quantum computing. 
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Figure Captions 



Fig. 1. The population difference s{t) (dashed hne) and phase difference x{t) (sohd hne) 
as functions of time, measured in units of a~^, for the zero detuning e = 0, initial conditions 
Sq — —1, Xq = 0, and varying amphtude of the resonant field: (a) b — 0.470; (b) b — 0.495; 
(c) b = 0.501; (d) b = 0.7; (e) 6=1. 

Fig. 2. Temporal behaviour of the population difference s{t) (dashed line) and phase 
difference x{t) (sohd line), under the same consitions as for Fig. 1, but for the negative 
parameter b <0 taking the values: (a) b = -0.40; (b) b = -0.49; (c) b = -0.6; (d) b = -1. 

Fig. 3. The population difference s{t) (dashed line) and phase difference x{t) (solid line) 
as functions of dimensionless time for £ = and the initial conditions sq = —0.8, xq — 0, 
with varying the parameter b as: (a) b — —1; (b) b = —0.85; (c) b — —0.8001 (slightly below 
the critical point bd); (d) b = —0.7999 (slightly above the critical point bd)', (c) b = —0.5; 
(f) b = 0.18; (g) b = 0.1999 (slightly below the critical point 602); (h) b = 0.2001 (slightly 
above the critical point 6c2); (i) b = 0.5; (j) 6=1. 

Fig. 4. Temporal behaviour of s{t) (dashed line) and x{t) (solid line) for the initial 
conditions Sq = —1, Xq = 0, fixed b = —1, and a nonzero detuning: (a) s = —0.05; (b) 
£ = -0.1. 

Fig. 5. Phase portrait on the plane {s.x} for a finite detuning e — 0.1 and different 
values of the pumping-field amplitude: (a) b = 0.51; (b) b = 0.8. 
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